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Abstract. Anisotropic flow is recognized as one of the main observables providing 
information on the early dynamics in heavy-ion collisions. The large elliptic flow 
observed at RHIC is considered to be evidence for almost perfect liquid behavior of 
the strongly coupled Quark Gluon Plasma produced in the collisions. In this report 
we review our current understanding of this new state of matter and investigate the 
predictions for anisotropic flow at the LHC. 
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1. Introduction 

Relativistic heavy- ion collisions are a unique tool to study the Equation of State (EoS) 
of extremely hot and dense matter under controlled conditions. The system of hot and 
dense matter created in a heavy-ion collision will expand and cool down. During this 
expansion the system probes a range of energy densities and temperatures, and possibly 
different phases. The collective dynamics of the system is revealed in the particle 
yield as function of transverse momentum, which is characterized by the temperature 
and transverse flow velocity of the system at freeze-out. Because heavy-ions are not 
point-like, the created system in non-central collisions has an azimuthal anisotropy 
in coordinate space which translates, due to multiple interactions, into an azimuthal 
anisotropy in momentum space. The elliptic flow t> 2 is defined as the second coefficient 
of a Fourier expansion of this azimuthal anisotropy pp. This elliptic flow is particularly 
sensitive to the early dynamics of the collision, for a recent review see [2]. 

Generally speaking, large values of elliptic flow are considered to be signs of 
hydrodynamic behavior as was first put forward by Ollitrault [3]. In ideal hydrodynamics 
t>2 is proportional to the spatial eccentricity with a magnitude which depends on the 
EoS. The spatial eccentricity is defined by 

= (y 2 - x 2 ) 

(y 2 + x 2 ) 
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Figure 1. (Color on-line) (a) Elliptic flow (solid points) as a function of centrality 
n ch/ n max- The open rectangles show a range of values expected for i> 2 in the ideal 
hydrodynamic limit [4: . (b) Elliptic flow of pions and protons as function of transverse 
momentum [5 . The curves are hydrodynamical model calculations using a hadron gas 
EoS (dashed curve) and an EoS which incorporates the QCD phase transition (full 
curve) . 



where x and y are the spatial coordinates of the colliding nucleons in the plane 
perpendicular to the collision axis and where the brackets denote an average. In practice 
e is not a measured quantity but obtained from Glauber or Color Glass Condensate 
(CGC) model calculations. The successful description of the measured elliptic flow [UG2] 
at RHIC by ideal hydrodynamical models, as shown in Fig. [TJ lead to the conclusion 
that the system created at RHIC behaves like a perfect liquid [6]. 

2. The perfect liquid? 

Also from Fig. [2] it is seen that the measured v 2 scaled by the spatial eccentricity versus 
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Figure 2. (Color on-line) Compilation of i^/e data [7J [S] versus particle density at 
midrapidity. Green lines indicate ideal hydrodynamic predictions for AGS, SPS and 
RHIC collisions energies 
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particle density reaches the expected ideal hydrodynamic values but only for the more 
central collisions at the highest RHIC center of mass energy E] . Discrepancies are 
observed for peripheral collisions, lower energies, and regions away from mid-rapidity 
which indicates that here the elliptic flow has significant non-ideal hydro contributions. 

Much of this discrepancy can be explained by incorporating the viscous contribution 
from the hadronic phase P, [101 El H2] • However, we expect that also the hot and dense 
phase must deviate from an ideal hydrodynamic description. Kovtun, Son and Starinets 
(KSS) [13], showed that conformal field theories with gravity duals have a ratio of shear 
viscosity 77 to entropy density s of, in natural units, rj/s = 1/4tt. They conjectured 
that this value is a lower bound for any relativistic thermal field theory. In addition, 
Teaney [14J pointed out that very small shear viscosities, of the magnitude of the bound, 
would already lead to a significant reduction in the predicted elliptic flow. 

Teaney [H] also pointed out that the reduction of elliptic flow due to non-zero rj/s 
would become larger at higher transverse momenta. In Fig. [3] ideal hydro calculations 
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Figure 3. (color online) Scaled elliptic flow t^/e versus transverse momentum for 
eccentricities used in ideal hydrodynamics (left) and CGC (right) compared to ideal 
hydrodynamics [T5J [TB] . 

are compared to STAR measurements of v^je as function of transverse momentum 
for different centralities. In the figure the initial eccentricity is calculated by either a 
Glauber model (left panel) or a CGC model [T7j (right panel). The figures show that 
with CGC initial eccentricities the magnitude olv%je is lower than with Glauber initial 
eccentricities [12]. 

The ideal hydro curves in Fig. [3] correspond to calculations [18J using an initial 
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eccentricity from a Glauber model. As was shown already in Fig. [2] the elliptic flow 
in the more central collisions reach the ideal hydrodynamical values, provided that 
Glauber initial conditions are used. It is seen in the right panel of Fig. [3] that this 
is not true anymore when CGC eccentricites are used which clearly illustrates that 
constraints on the initial conditions are crucial. Nevertheless, for both Glauber and 
CGC initial eccentricities, it is seen that going from peripheral to central collisions 
the ratio Vije increases and that the transverse momentum where v% is maximal rises. 
Both these observations are consistent with a decrease of r\js from peripheral to central 
collisions [H]. 

One approach to quantify the contribution from rj/s is to incorporate viscous 
corrections into the hydrodynamic model. A drawback of this approach is that solving 
viscous relativistic hydrodynamical equations is not straight forward and that the results 
are sensitive to not only rj/s but as in ideal hydrodynamics also to, among others, the 
initial spatial eccentricity [12] and the EoS [19] . Nevertheless, tremendous progress has 
recently been made in this area (2QI EH [221 EH ESI EHl EH ESI [29l [30]. 
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Figure 4. (Color on-line) (a) Elliptic flow versus transverse momentum using 
CGC initial conditions in hydrodynamics with viscous corrections. The STAR 
measurements (open circles) have been approximately corrected for fluctuations and 
nonflow. (b) Fit to v 2 /e versus particle density in terms of Knudsen number for Cu+Cu 
and Au+Au [17]. 

Figure [4^i shows, as an example, the measured transverse momentum dependence of 
the elliptic flow for charged particles compared to a viscous hydro model calculation |26j . 
Here, the initial eccentricity is based on a CGC model calculation and results are shown 
for various values of r]/s. It is again clear from Fig. [4] that the difference between ideal 
hydrodynamics (full curve) and the measured flow is substantial and that this difference 
increases with transverse momentum. Good agreement with the data is obtained when 
rj/s is about twice the KSS bound of 0.08. 

Motivated by transport calculations, another promising approach to quantify the 
possible discrepancy with ideal hydrodynamics is a description of v-ije versus particle 
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K = ac s 



density in terms of the Knudsen number (TTJ 151] . This Knudsen number is the mean 
free path of the constituents divided by the system size, and is defined by [T7J [31]: 

i dNy 1 

' S dy ) 

Here a is the constituent transport cross section, c s the velocity of sound and S the 
transverse area Ait^J (x 2 ) (y 2 )- 

The scaled elliptic flow is parameterized as [T7J [31] : 

V -l = h (1) 

e 1 + K/K ' {) 

where h is the ideal hydrodynamic limit of V2/S and Kq = 0.7 is a constant determined 
from transport calculations [TTJ . This parameterization yields correctly the hydro limit 
for K — > and the low density limit [32] for K — > 00: 
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Figure |4p shows i^/e measured by PHOBOS in Cu+Cu and Au+Au collisions at 
\/ s iviv = 200 GeV, using Glauber initial conditions for e. The curve in the figure shows 
a fit of Eq. [T] to the data with h and a as free parameters. Input to the fit are the 
quantities dN/dy, e and S obtained from measurement and model calculations. Clearly 
the parameterization describes the data well over the whole range of particle densities 
and shows no sign of saturation. Therefore the authors of the fit conclude in [17] that 
even for the most central collisions at the highest RHIC energy the deviations from 
ideal hydrodynamics are as large as 30%. The magnitude of r\j s is obtained from the 
fitted value of o combined with an additional estimate of the temperature where 
develops. The value of rj/s is found to be a few times the KSS bound, which is roughly 
in agreement with estimates from viscous hydrodynamical model calculations. 



3. Collective flow at the LHC 



Most theoretical models [33] predict an increase of v 2 at the LHC compared to RHIC. 
In the following we will describe estimates based on (ideal) hydrodynamics, the low 
density limit and the Knudsen parameterization. 

Figure [5] shows predictions of t^/e versus center of mass energy from a model 
that combines ideal hydrodynamics and a hadron cascade [31]. The lower curve (full 
triangles) shows the elliptic flow from the partonic phase, assuming ideal hydrodynamics. 
The upper curve (full squares) shows this elliptic flow with the contribution from the 
hadronic phase added, here the hadronic contribution also is calculated using ideal 
hydrodynamics. When using a hadron cascade instead of hydrodynamics to describe 
the evolution of the hadronic phase the predictions are obtained as shown by the full 
circles in Fig. [5j It turns out that this combination gives the best description of the 
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Figure 5. (color online) v-i versus center of mass energy from a ideal hydro + hadron 
cascade model [3J]- 



data (not shown). It is clear from these calculations that the contribution from the 
partonic phase is considerably larger at the LHC than at RHIC energies. It is also seen 
that the prediction from a full ideal hydrodynamic evolution yields the same answer as 
a hydro + hadron cascade evolution at LHC energies. This indicates that the viscous 
contribution from the hadronic phase is relatively unimportant at the LHC. 




Figure 6. (color online) (a) V2 versus 77 measured by PHOBOS at RHIC, using 
longitudinal scaling to make a prediction for the LHC [36]. (b) Elliptic flow t>2 as 
function of centrality N part at mid-rapidity for Pb-Pb collisions at the LHC 37J. 

In the low density limit elliptic flow is proportional to the particle multiplicity. It 
has been shown by PHOBOS that particle multiplicity as function of center of mass 
energy is well described using longitudinal scaling [35J. It is therefore expected, in 
the low density limit, that also t>2 exhibits longitudinal scaling, provided that the 
eccentricity is kept constant. The longitudinal scaling is demonstrated in Fig. [6^1 where 
measurements of t>2 at various RHIC center of mass energies are plotted versus the 
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pseudorapity rj at the LHC. Linear extrapolation to mid-rapidity yields an increase of 
about 50% in the magnitude of v 2 compared to RHIC. This increase in v 2 is larger than 
predicted by most other models. 

The third prediction presented here is based on the parameterization in terms of 
the Knudsen number (Eq. fl|. The two lower curves in Fig. [6|d show the centrality 
dependence of v 2 assuming either a fixed coupling or running coupling evolution of the 
saturation scale Q s . The corresponding transport cross sections are a = 5.5 and 3.3 mb 
respectively. The first value corresponds to RHIC assuming CGC initial conditions. The 
upper two curves are prediction based on "ideal hydrodynamics" (assuming v 2 /s = h). 
For fixed coupling it is seen that the Knudsen parameterization closely approaches the 
ideal hydro curves except for the most peripheral collisions. The ideal hydrodynamical 
limit h is a free parameter therefore the predicted magnitude of v 2 itself is not fixed 
(however all four curves move together changing h) . For the figure the previously found 
RHIC value of h = 0.22 is used. The centrality dependence of v 2 does not depend on h 
and, as can be seen in the figure, the measurement of this centrality dependence at the 
LHC will allow us to distinguish between ideal hydro and the Knudsen parameterization. 

4. Summary 

At RHIC the observed large elliptic flow provides compelling evidence for strongly 
interacting matter which appears to behave like an almost perfect fluid. To quantify 
this the significant effects of the viscous corrections need to be calculated. At y/s NN = 
200 GeV, based on different model assumptions, calculations oirj/s have been performed 
and are found to be in approximate agreement, ranging for the more central collisions 
from 1-4 times the KSS bound. Even though this is encouraging, it is important to 
realize that these results are all obtained based on a small set of initial conditions. The 
estimated range of rj/s values would certainly be affected if the initial conditions were 
to be sufficiently different, e.g. due to strong initial flow fields. 

To quantify rj/s of the partonic fluid requires knowledge of the relative contributions 
from the partonic and hadronic phase. Detailed comparison of identified particle elliptic 
flow, in particular for particles which have different hadronic cross sections, should allow 
us to constrain this experimentally. In addition, comparing the center of mass energy 
dependence of v 2 in which the relative contribution from the partonic and hadronic 
phase to v 2 varies should provide additional constraints. Elliptic flow measurements at 
the LHC will provide an important contribution to this and in addition a decisive test 
which of the currently successful model descriptions is more appropriate. 
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